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We consider the early Universe scenario which allows for production of non-Gaussian 
curvature perturbations at small scales. We study the peculiarities of a formation of 
primordial black holes (PBHs) connected with the non-Gaussianity. In particular, we 
show that PBH constraints on the values of curvature perturbation power spectrum 
amplitude are strongly dependent on the shape of perturbations and can significantly 
(by two orders of magnitude) deviate from the usual Gaussian limit Vr < 10~ 2 . We give 
examples of PBH mass spectra calculations and PBH constraints for the particular case 
of the curvaton model. 
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1. Introduction 

As is well known, in models of slow-roll inflation with one scalar field the curvature 
perturbation originates from the vacuum fluctuations during inflationary expansion, 
and these fluctuations lead to practically Gaussian classical curvature perturbations 
with an almost flat power spectrum. However, it is well known also that both these 
features are not generic in the case of inflationary models with two (or more) scalar 
fields: such models can easily predict adiabatic perturbations with, e.g., a "blue" 
spectrum and these perturbations can be non-Gaussian.El 

Possibilities for appearing of non-Gaussian fluctuations in inflationary models 
with multiple scalar fields had been discussed long agoPEl The time evolution of 
the curvature perturbation on superhorizon scales (which is allowed in multiple- 
field scenarios^) implies that, in principle, a rather large non-Gaussian signal can 
be generated during inflation. According to the observational dataj^the primordial 
curvature perturbation is Gaussian with an almost scale-independent power spec- 
trum. So far there is only a weak indication of possible primordial non-Gaussianity 
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[at (2 — 3)tr level] from the cosmic microwave background (CMB) temperature in- 
formation data (see, e.g., Ref. |7|). However, non-Gaussianity is expected to become 
an important probe of both the early and the late Universe in the coming yearsP 
The second important feature of predictions of two-field models is that these 
models can lead to primordial curvature perturbations with blue spectrum (for 
scales which are smaller than cosmological ones) and, correspondingly, can predict 
the primordial black hole (PBH) production at some time after inflation. In this case, 
PBHs become a probe for the non-Gaussianity of cosmological perturbationsP^l 
The results of PBH searches can be used to constrain the ranges of early Universe 
model parameters. 

There are several types of two-field inflation scenarios in which detectable non- 
Gaussianity of the curvature perturbation can be generated: curvaton modelsP^^ 
models with a non-inflaton field causing inhomogeneous reheating] curvaton- 
type models of preheating (see, e.g., Ref. HU and references therein), models of 
waterfall transition that ends the hybrid inflationP3tH| In these two-field models 
the primordial curvature perturbation £ has two components: a contribution of 
the inflaton (almost Gaussian) and a contribution of the extra field. This second 
component is parameterized by the following waj^ 

C x (x)=a X (x)+ X 2 (x)-( X 2 ). (1) 

If a = 0, one has a x 2 -model. Obviously, the quadratic term can't dominate in £ on 
cosmological scales where CMB data are available. It can, however, be important 
on smaller scales. 

In the present work we study the predictions of the PBH production and cor- 
responding PBH constraints for the curvaton model. The potentially large non- 
Gaussianity in this model is connected with the fact that the predicted magnitude 
of the curvature perturbation is proportional to a square of the non-inflaton (cur- 
vaton) field. The blue spectrum in the curvaton model is due to, e.g., supergravity 
effects leading to the large effective mass of the curvaton.^ 

The main attention in the present paper is paid to a study of probability dis- 
tribution function (PDF) of the curvature perturbation and the shape of the black 
hole mass function, with taking into account of the non-Gaussianity. The first gen- 
eral study of PDF of the curvature perturbation in curvaton model was carried out 
in Ref. [23 

PBH production in curvaton scenario was studied in recent works 26, 27 (without 
considering the non-Gaussian effects). The approximate PBH constraints on the 
curvature perturbation power spectrum in the curvaton model were obtained in 
Ref. EH 

The plan of the paper is as follows. In Sec. [2] we study the process of PBH 
production in the case when the primordial curvature perturbations are strongly 
non-Gaussian. We calculate the PDF function for the % 2 -model and use it for a 
calculation of PBH mass spectrum via Press-Schechter mechanism. In Sec. [3] we 
discuss the possible production of PBHs in the curvaton model and the correspond- 
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ing cosmological constraints that can be obtained. Our conclusions are given in Sec. 

SI 

2. PBH mass spectrum in the case of non-Gaussian curvature 
perturbations 

2.1. PDFs for x 2 '-model 

Generally, in x 2 -model the connection between curvature perturbation £ and the 
square of the (non-inflaton) scalar field perturbation value y 2 is 

C = ±A( X 2 < x 2 )), A > 0. (2) 

The distribution of \ is assumed to be Gaussian, i.e., 

p x ( X ) = ^=e"^, a 2 ^( X 2 ). (3) 



The form of the PDF depends on the sign in front of the right-hand side of Eq. ([2]) . 
If the sign is negative, i.e., 

C = -A( X 2 -(X 2 )), (4) 
it is convenient to introduce the notation 

Cmax = A(x ), C — Cmax, (5) 

and the PDF is given by the formula^ 

PdO =Px 



e^c^, (6) 



which is just a ^-distribution with one degree of freedom, with an opposite sign of 
the argument, shifted to a value of C, max . 

In a case of the positive sign in Eq. one has, correspondingly, 

C = A( X 2 -(x 2 )), (7) 

and, introducing the notation 

( mm = ~A( x 2 ), C>C™„, (8) 
one obtains for the PDF the expression 

PdO= . ) == e^tt. (9) 

The case when the sign is negative is realized in the model of hybrid inflation 
waterfalP^"^ and was studied in detail in Ref. [29] In this case, x ls the perturbation 
of the waterfall field. 
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In our present case, when x is the perturbation of the curvaton field, the sign is 
positive. The variance of the PDF of the £ field distribution is 

oo 

(C 2 )= J ( 2 pdC)d( = 2d„. (10) 

Cmin 

This variance is connected with the curvature perturbation spectrum "P^ through 
the expression 

(C 2 )=« 2 = ■ (ii) 

The distribution function ^ can be written in the form 

PC(0 = —P (-) = —Piy\ P(") = -=L=e-^+^\ (12) 

Here, the ratio v = C/°"C is introduced. The first central moments of the PDF of 
the £ field are given by 

(0=0; (C 2 )=2A 2 ( X 2 ) 2 ; (C 3 ) = 8A 3 ( X 2 ) 3 ; (C 4 ) = 60A 4 ( X 2 ) 4 , (13) 

where (x 2 ) is a variance of the X -field power spectrum, 

(X 2 ) = / P x (*)f. (14) 

We will use also the central moments of p{y) distribution ("reduced central mo- 
ments"). They are given by the expression (( n )/a^. In particular, for the skewness 
and kurtosis one has, respectively, 

S -<£>--<£>- (15) 

a 3 (£2)3/2' ^ (£2)4/2- ^ 

For the following, we need the expression for the PDF of the smoothed curvature 
fluctuations, i.e., instead of Eq. ([7]) we must use the smoothed C field, 

Ch(x) = A J dyW>(|x - y|/i?)x 2 (y) - A( X 2 ) J dyW(\x - y\/R), (16) 

where W is the window function. We will use the Gaussian form of this function, 
in this case its Fourier transform is W(kR) = exp(— k 2 R 2 /2). The expressions for 
the central moments of the corresponding PDF, P(;,r, had been derived in Ref. |321 
The second and third central moments of pq,r are given by the formulas 

(Cfl) - -JrgA 3 J dkdk'P x (k)P x (k')W 2 (\k + k'\R), (17) 

(C « } = (2^/ dkdk ' dk " P x(k)P x (k')P x (k") x 

x W(|k + k'\R)W(\k' - k"\R)W(\k + k"\R), (18) 
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where P x (k) is the power spectrum of the x held, 

p x (k) = 2 -^-r x (k). (19) 

Now, we suppose, that the PDF of the smoothed ( field can be presented in the 
factorized form [as in Eq. (|12|)]. 

P<,R = —T^P ( —Td\ ) = —7n\P( v R)i ( 20 ) 

^)HC£> 1/2 - (21) 

If Eq. (|20[) is approximately correct, the central moments of p(vr) are weakly de- 
pendent on the smoothing scale R. It had been shown in Refs. |33j |34j |35j [36] that 
it is really so, for a wide range of scales. In particular, SetcS^ showed that, if the 
spectrum V x has a power form, V x ~ fc* x , the scale dependences of the variance 
and third central moment of the £r field are: 

(Cl)~ir 2 \ (d)~ir 3 **. (22) 

Correspondingly, the scale dependence is canceled in the expression for the skewness 
parameter, Sr = ((r) / ((r) 3 ^ 2 ■ Moreover, it appear^ (and it is most essential for 
our case) that, if t x < 1, the value of Sr is quantitatively close to S [which is equal 
to V&, as follows from Eq. (fl"3|) ]. The analogous check had been performed in Ref. 
[M]for the next central moment (kurtosis). 

It follows from this analysis that the shape of vr distribution is close to the 
shape of v distribution, and the smoothing effects enter only through the value of 
the variance, {(r) 1 ^ 2 = &C,{R)- 

As we will see in the next Section, in our case the spectrum of the x-ficld has a 
power form, V x ~ k tx , and t x is of order of 1. Therefore, basing on the above cited 
works, we will use Eq. (|20|) for the PDF of the smoothed ( field, with p(vr) having 
the same form as p{v). 

The form of the distribution ^ for cr^ = 2 x 10~ 4 is shown in Fig. [TJ It is seen 
that in this particular case the probability to reach £ ~ ( c ~ 1 (which is, as is well 
known, required for the PBH formation) is ~ 10~ 20 or so, i.e., roughly of the same 
order of magnitude as PBH constraints on energy density fraction of the Universe 
contained in PBHs at the time of their formation, Ppbh(Mbh), in PBH mass range 
M BH - (10 10 - 10 20 ) gPlt follows from this Figure that the value of V ( (k) - 10~ 4 
is already enough for producing an observable amount of PBHs in this model (this 
is in agreement with the estimates of Refs. l30l [28| . 

2.2. PBH mass spectrum in Press-Schechter formalism 

For a derivation of PBH mass spectrum and PBH constraints we will use the Press- 
Schechter formalism!^ We will follow the approach of Refs.[28l[29] and more recent 
works 39, 40 working directly with the curvature perturbation rather than with the 
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Fig. 1. The form of the distribution for = 2 X 10 4 . Dashed lines show the considered 
values of £ c (0.75 and 1, see text). 

density contrast. In the Press-Schechter formalism, the energy density fraction of 
the Universe contained in collapsed objects of initial mass larger than M is given 

by 



— / Mn(M)dM 
Pi J 

M 



p c (()d( = P((>t c ;R(M),t i ), 



(23) 



where function P in right-hand side is the probability that in the region of comoving 
size R the smoothed value of £ will be larger than the PBH formation threshold 
value, n{M) is the mass spectrum of the collapsed objects, and pi is the initial 
energy density. The parameter £ c in Eq. (|23|) is the threshold of PBH formation 
in the radiation-dominated epoch, which is to be taken from gravitational collapse 
model. For estimates, in the following we will use two values: £ c = 0.75 and £ c = 1 
(corresponding to the PBH formation criterion in the radiation-dominated epoch: 
5 > S c , with 6 C = 1/3 and 8 C = 0.45, respectively).^ 

The horizon mass corresponding to the time when fluctuation with initial mass 
M crosses horizon is (see Ref. |4T|) 

M h = MV 3 M 2/3 , 



where Mi is the horizon mass at the moment t 

47T 



Mi 



tlp l 



(24) 



(25) 



and U is the time of the start of the radiation era, pi is the energy density at this 
time. If the reheating is short, U coincides with the time of the end of inflation. 
In this case, there is a connection between Mi and the reheating temperature (see, 
e.g., Ref. EI]): 

3 



Mi w 0.038- 



' pi 



fJ 



l/2r 



100. 



(26) 



RH 
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For simplicity, we will use the approximation that mass of the produced black 
hole is proportional to horizon mass, namely, 

Mbh = f h M h = f h Mj /3 M 2 / 3 , (27) 

where f h (1/3) 1 / 2 = const. 

Using (1231) and (|2?p , for the PBH number density (mass spectrum) one obtains^ 

dP 



ubhIMbh) = — 



4vrV 1/3 dP f h Pi /3 M} /a 



dR 



ai M 2 BH 



(28) 



where di is the scale factor at the moment U. The derivative dP/dR (where P is 
the function defined in Eq. (|23p) is given by the expression 

dP Cc da c (R) 

dR-a c (R) dR PC ' fl(Cc) (29) 

(in derivation of this equation the PDF form of Eq. (|20[) was used). 

The dependence of the PBH mass spectrum on the curvature perturbation spec- 
trum enters just through the factor dP /dR. 

Introducing the PBH formation time, t — t e (see Sec. [3]), we can calculate the 
energy density fraction of the Universe contained in PBHs, at the moment t e (at 
this moment the horizon mass is equal to Mh(t e ) = M[)W^ 

i ( m s \ 1/2 f 

n PBH (Mf) « - I -^A J J n BH {M BH )M 2 BH d\nM BH « 




n B H{M BH )M% H 



PiM,, 



1/2 



n BH (M BH ) \ MnH=M rnin ■ (30) 



Here, M B % H ^ is the minimum mass of the PBH mass spectrum, Mg]ff = fhM[ . It is 
well known that for an almost monochromatic PBH mass spectrum £Ipbh(m[) 
coincides with the traditionally used parameter Ppbh- Although all PBHs do 
not form at the same moment of time, it is convenient to use the combina- 
tion M i p^ 1 Mg£n B ij(M B H) to have a feeling of how many PBHs (with mass 
~ Mbh) actually form, i.e., to use the estimate following from (130|) : 

Mr 1 !" 1 p^M^lnBHiMBH) w Ppbh ■ (31) 



3. PBH constraints in the curvaton model 

Curvaton is an additional to the infiaton scalar field that can be responsible (partly 
or fully) for generation of primordial curvature perturbationsP^^H^ This field can 
also be a source of PBHs, as discussed in Refs. [42J l30l 

In this work, we only consider the case of a strong positive tilt (the possibility 
discussed in Refs. [TJ II 5|) of the curvaton-generated perturbation power spectrum. 
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At the same time, it is assumed that inflaton is responsible for generation of per- 
turbations on cosmological scales (see Fig. [2] for an illustration) . 

The curvaton held generates cosmological perturbations in two stagesPHK^ 

(i) Quantum fluctuations of the curvaton during inflation (at time of horizon 
exit) become classical, super-horizon perturbations. 

(ii) In the radiation-dominated stage, the curvaton starts to oscillate (this hap- 
pens at the time when Hubble parameter becomes of order of curvaton's effective 
mass, H ~ to). The Universe at this stage becomes a mixture of radiation and mat- 
ter (the curvaton behaves as a non-relativistic matter in this regime) . The pressure 
perturbation of this mixture is non-adiabatic and the curvature perturbation is thus 
generated. One obtains, approximately, the expression (see, e.g., Ref. 

' / ' \ 2 

rt7 osc r„ i r ( °osc 1 c 2 



C(t,x) = -2S^ + - -SSi Sai, (32) 

la osc 4 \<7osc/ 

where r is the density parameter, r = 4p (T /(4p r + 3p ff ) (p r is the energy density 
of radiation after inflation), er osc is the value of the curvaton field at the onset of 
oscillations. The initial value for the curvaton field, cr*, is set by inflation. The 
derivative in Eq. (|32[) is taken with respect to the field value during inflation, ov 
The term containing the second derivative, <r" s[ , is neglected. It is assumed that 
r < 1. 

Assuming zero average value for the curvaton field (i.e., working with the max- 
imal boxP^), we keep in Eq. (|3"2")) only the second term, 



r 



2 



C(t,x) = -I^c So 2 (33) 

4 \ Cr osc J 

In this case we have the % 2 -model (because the perturbations <5cr* are assumed to be 




Fig. 2. A sketch that illustrates a relation between curvaton-generated and inflaton-generated 
curvature perturbation power spectra for the scenario of PBH production that we consider. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
log 10 (M BH /M BH min ) 



Fig. 3. Examples of PBH mass spectra calculated for the curvaton model, for several sets of 
parameters n and £ c . For all curves, P' = 4x 10 -4 . The mass M^ff = fhMh(te) and t e is the 
moment of time when perturbation with comoving wave number k e enters horizon. 



Gaussian). In notations of Sec.[2j one has A = j (f 22 ^ , — X- The fluctuations 
are strongly non-Gaussian which is not forbidden on small scales. 

The curvaton-generated curvature perturbation spectrum can be writterP^^H 
as (using the Bunch-Davies probability distribution for the perturbations of the 
curvaton fielcPEsl) 



V, 



1/2 



2 

3^ 



1 



2^ 



(34) 



^ 0"osc J \to 

where fl a — p a / P ~ T is the relative curvaton energy density at the time of its decay 
(it depends on the curvaton decay rate) and 

2ml 



(35) 



is the tilt of the curvaton field spectrum, V a ~ k ta , m* is the effective mass of the 
curvaton field and is the Hubble parameter during inflation. 

It is rather natural (see Ref. [1] and, e.g., Ref. [47j which considers the models 
of chaotic inflation in supergravity) that t a ~ 2/3 which corresponds to a blue 
perturbation spectrum with the spectral index 

n = 1 + 2t a - 7/3 (36) 

(such a situation is shown in Fig. [2]). For the following, we parameterize the spectrum 
in a simple form 

-l 



(37) 



and will treat 7-^, n, and k e as free parameters. Note that k e does not, in general, 
coincide with the comoving wave number corresponding to the end of inflation. 
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logio M h [g] 

30 28 26 24 22 20 18 16 14 12 10 




8 10 12 14 16 18 

log 10 k e [MpcT 1 ] 



Fig. 4. The limits on the maximum value of curvature perturbation power spectrum from 
PBH non-observation, for the curvaton model, for several sets of parameters n and £ c . The for- 
bidden regions are above of the corresponding curves. Mh is the horizon mass at the moment 
of time when the perturbation with comoving wave number k e enters horizon. The spectrum is 
parameterized as \37). 



Rather, it is the scale entering horizon at the time when £ is created^ (we assume 
that it is created instantaneously due to the fast curvaton decay). 

Using the Eq. (|28l) , we can calculate PBH mass distributions that are generated 
for a particular set of parameters (n, V®, etc.) and then compare the resulting 
Ppbh (Eq. ([51)) with the known limits (from, e.g., Ref . [37J) . The example of PBH 
mass spectrum calculation is shown in Fig. [31 It is seen from this Figure that PBH 
abundances strongly depend on the particular choice of £ c . 

In the calculation of PBH mass spectra we took into account the fact that PBHs 
start to form only after the moment of time t = t e when curvaton decays and £ is 
created,^ and the scale k e enters horizon just at the same moment t e . Thus, the 
minimal PBH mass in our model, taking into account Eq. (|27p , is 



MTh = hM h (t e ). (38) 

This minimal mass corresponds to the vertical line in Fig. [3J 

The resulting constraints on parameter V® (for two values of the spectral index 
n) following from datcP^ on PBH non-observation are shown in Fig. [4] Here, the 
connection between k e and Mh is given by (see, e.g., Ref. |48|) 

^ 2 x 10 23 ! 

The obtained constraints can also be transformed into limits on particular curvaton 
model's parameters, such as f2 CT . For example, comparing Fig. [4] and Eq. (|34[) . one 
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obtains, roughly, 

ft* « (7> c °) 1/2 < 10~ 2 , (40) 

which is already a useful constraint. It is a subject of our further study to get more 
exact limits for particular parameter sets of the model. 

4. Conclusions 

Primordial black holes can be used to probe perturbations in our Universe at very 
small scales, as well as to study other problems of physics of early stages of the 
cosmological evolution. We have considered the PBH formation from primordial 
curvature perturbations produced in the curvaton model. This model predicts the 
production of strongly non-Gaussian perturbations, and non-Gaussianity was taken 
into account in the calculation of PBH mass spectrum (in Press-Schechter formal- 
ism). Limits on the values of perturbation power spectrum as well as approximate 
constraints on inflation model parameters were obtained. The constraints on the 
curvature perturbation spectrum amplitude follow from Fig. 2J It had been shown 
in our previous worlP^ that the constraint on Vq in the case when PDF of the 
curvature perturbation is given by Eq. (j4)) (this case is realized, e.g., in the hybrid 
waterfall model 2 ^) is very weak, "P^ < 1. In contrast with this, the corresponding 
constraint in the curvaton model is much stronger, Vq < (10 -4 — 10~ 25 ), depend- 
ing on the value of the spectral index and PBH mass. It is important to note that 
the latter constraints are more strong than those following from one-field inflation 
models, in which < 10 -2 (see, e.g., Ref. |48|) . The PBH constraints obtained in 
this work confirm the estimates given in Refs. 1301 l2"51 
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